Journal of Modern Optics 

Vol. 00, No. 00, DD Month 200x, 1-11 



1^ Compression and localization of an atomic cloud in a time dependent optical lattice 
O . 

Q . Will Williamst & M. Saffman 

Department of Physics, 1150 University Avenue, University of Wisconsin, Madison, WI, 53706 

(Received 00 Month 200x; in final form 00 Month 200x) 



(N 



> 
O 

We analyze a method of compressing a cloud of cold atoms by dynamic control of a far off resonance optical lattice. We show that by 
reducing the lattice spacing either continuously or in discrete steps while cooling the atoms with optical molasses large compression factors 
can be achieved. Particle motion in the time dependent lattice is studied numerically using a three dimensional semiclassical model. Two 
experimentally realistic models are analyzed. In the first we continuously vary the lattice beam angles to compress atoms initially in a 
^ Gaussian distributed cloud with standard deviation of 250 /im into a single site of a two-dimensional lattice of area A ~ 35 x 35 A^, 

« with A the wavelength of the lattice beams. This results in an optical depth for an on-resonant probe beam > 80 which is an increase 
^ * by a factor of about 1800 compared to the uncompressed cloud. In the second approach we use a discrete set of lattice beam angles to 

I decrease the spatial scale of the cloud by a factor of 500, and localize a few atoms to a single lattice site with an area A ~ A^. 
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■ 1 Introduction 

During the last few decades there has been intense development of methods to cool and localize atomic 
^ samples with optical fields. Conservative trapping potentials can be created using far off resonance light 
^ ^ beams that are focused to a small spot [1,2]. These far off resonance traps have become widely used as 
a starting point for achieving Bose Einstein condensation [3-5], as well as for experiments with single 
atoms [6-8]. There is currently much interest in techniques for loading a large number of atoms at high 
^ densities into optical traps [9] and in methods of controlling the loading of an optical lattice [10]. These 
^ capabilities are important for creating large condensates, as well as highly localized atomic ensembles for 
^ quantum logic and communications [11]. The phase space density of atomic clouds can be increased using 
^ laser or evaporative cooling as well as by novel approaches such as asymmetrical optical barriers [12]. 
O Heating due to photon rescattering limits the phase space density that can be reached by optical cooling 
although spatially anisotropic samples [13] as well as optimized cooling techniques in optical lattices [14,15] 
have led to improved performance. Dynamically changing the optical potential in conjunction with optical 
cooling can also increase phase space density as in experiments which compressed an atomic cloud using 
O a trap created by rapidly rotating a single blue detuned beam [16, 17]. 
^ In this paper we use numerical simulations to study the performance of a novel protocol for transferring 
^> a three-dimensional distribution of atoms into a single well of a planar lattice. The anisotropic well is 
elongated perpendicular to the plane of the lattice which minimizes heating and repulsion due to photon 
rescattering [13] and is well suited for applications that require mode matching between a Gaussian laser 
. ^ beam and a pencil shaped atomic ensemble [11]. The technique can also be extended to a three dimensional 
^ lattice to create a spherically shaped localized cloud of atoms. 

The proposed method starts with a large period optical lattice which confines all atoms in a single 
trapping site. The angles between the beams are then increased to create wavelength scale lattice sites. 
We show that if the lattice beam angles are changed slowly, and the atoms are simultaneously cooled by 
optical molasses, it is possible to keep a large fraction of the atoms in a single highly localized trapping 
site. Atoms of temperature T are localized within a distance 6r w ^/UbT/U from the origin in a well of 
depth U that is created by a focused laser beam with Gaussian waist (radius to the 1/e^ intensity point) w. 
The atoms have a vibrational oscillation frequency oi u = {2/w)^/uJm, where m is the mass of an atom. 
Provided the well size is reduced on a time scale that is long compared to the oscillation frequency the 
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Figure 1. Arrangements for dynamic compression by varying the angles between lattice beams: a) The beam crossing angles are 
changed continuously using four pairs of galvanometer mounted mirrors, b) The beam crossing angles are changed in discrete steps 
using two SLMs. The additional cooling beams provide optical molasses during the compression. 



atoms will tend to be trapped in the center of the well and get compressed. Under this adiabatic condition 
optical molasses is not necessary to spatially compress the atoms but there will be no increase in phase 
space density. By combining dynamic lattice reshaping with optical cooling the atoms can be compressed 
to higher phase space density. 

In this paper we numerically explore two methods for manipulating atomic samples using dynamic 
optical lattices. The first method assumes a continuously variable angle between lattice beams, as in Fig 
[1^, to compress a large cloud in two dimensions to provide a large optical depth along the third dimension. 
Continuously variable angles could be obtained using, for example, galvanometer mounted mirrors or 
acousto-optic modulators. The second method uses a discrete set of beam angles, as could be implemented 
with one-dimensional spatial light modulators (SLM's) with N pixels, see Fig [lb. In the first approach 
we use relatively small compression factors of about 15 to create pencil shaped atom clouds with a large 
optical depth. In the second approach we can achieve a compression factor of about 500 and localize atoms 
to a single lattice site with area A ^ A^, with A the wavelength of the lattice beams. 

The rest of the paper is organized as follows. The theoretical and numerical model for the dynamical 
optical lattice are detailed in Sec. El Section [3] presents results of numerical simulations that validate the 
proposed compression technique for continuous variation of the beam angles, as well as discrete steps. A 
comparison between traveling wave and standing wave lattice geometries is also given. We conclude with 
Sec. m which gives a discussion of applications and limitations of the technique. 



2 Theory 

We describe a method to compress atoms in a single site of a two-dimensional lattice by dynamically 
changing the angles between the lattice beams. The conservative optical potential for a two-level atom 
including the effects of saturation is [18] 



1 + 



/(r) 



I, (1 + 4A2/72) 



(1) 



where /(r) is the optical intensity at position r, Ig is the saturation intensity, 7 is the natural linewidth, 
A = — is the detuning from resonance, uj is the optical frequency, and Ua is the atomic transition 
frequency. 

Consider two equal intensity plane traveling waves with wavevectors kj = k{cos0jX + sinO jy) where 
k = 27r/A, A is the wavelength, and 0j is the angle each beam makes with the x— axis. When 9i—-92—0 
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Figure 2. Left: Density plot of the interference pattern created by four traveling waves, as described in the text, with 6 = tt/S. Regions 
of higher intensity are shown in white. Right: The intensity along the y— axis as a function of 0. When ^ = 0, there is a plane traveling 
wave with equal intensity over all space. As increases, a periodic structure forms and the well area decreases. 

the angle between the beams is 29 and the intensity is 

= 4/oCOs2[7ry/A(^)] (2) 

where Iq is the intensity of each wave and the spatial period is A(^) = Js^jo)' achieve compression 
in one dimension 9 is increased from to 7r/2 radians. Doing so decreases the periodicity from infinity 
(a plane wave with uniform intensity) to ^, the periodicity of a single standing wave. To accomplish two 
dimensional compression another set of traveling waves are added with their bisector orthogonal to the 
bisector of the first pair as in the configurations shown in Fig. [H By arranging for the second pair of waves 
to be polarized orthogonally to the first pair we obtaiij^ 

I{x,y) = 4/0 {cos^ [7Tx/A{9)] + cos^ [7Ty/A{9)]} . (3) 

This intensity produces a two dimensional lattice pattern with periodicity of A in both the x and y 
directions. Each site of the lattice has an area of A{9) ^ = 4sin^(0) ' ^ increases from 0, the area 
of each well decreases as shown in Fig. [2l If the initial value of 9 is small enough the central well will be 
large enough to contain essentially all of the atoms. Any atoms that start in the central well will stay there 
during the compression of the wells as long as they are not energetic enough to travel over the intensity 
barrier into the next well. If some of the atoms escape the central well, the method can be repeated to 
obtain better compression. A similar protocol is also possible using standing waves as will be discussed in 
Section [3l 

The compression protocol thus involves the following sequence. We start with a cloud of cold atoms 
and superimpose a two-dimensional lattice of optical trapping beams created by two pairs of crossing 
beams propagating along x and y as shown in Fig. [H The initial value of = ^min is chosen small 
enough such that almost all the atoms are confined in the central lattice well. We then increase 9^ either 
continuously or in discrete steps, which spatially compresses the atoms but also heats them since the 
well area decreases. The increase in temperature is removed by applying optical molasses during the 
compression, either continuously, or for a time Tc between each step. The compression sequence is applied 
until the atoms have been localized into a lattice site with small area created by the maximum crossing 
angle 9 6>max. 




■"^ Superposition of the orthogonally polarized pairs of beams creates a field that has a spatially varying polarization state. In a real 
multilevel atom the optically induced potential depends on the polarization when there is a nonzero vector or tensor polarizability. We 
assume in this paper that polarization effects are negligible so that the potential is accurately described using the intensity of Eq. (|3)). 
This is true for alkali atom ground states provided the wavelength is sufficiently far detuned from the strong D lines. 
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3 Numerical Results 

Atomic compression was simulated numerically using Newtonian dynamics. The equation of motion mr = 
—VU was used to update the velocity and position of each atom during the compression protocol. Atomic 
collisions were not accounted for. Cooling by three dimensional molasses was simulated by testing the 
probability that a photon is scattered at every time step using an approximation to the three dimensional 
distribution of the optical field as described in [19]. If the atom scattered a photon, the atom's momentum 
was changed by hk along the direction of the appropriate cooling beam for absorption and then in a random 
direction in three dimensions for emission. The dipole emission pattern of a real atom was approximated 
by a spherically uniform emission probability. The atomic parameters were chosen to correspond to the 
D2 transition of Cs (6^51/2 — at A = 852 nm, with excited state decay rate 7 = 27r x 5.22 MHz 

(this was reduced by a factor of 3 as explained below). The intensity of each cooling beam was set to 
I /Is = 0.1 (1.09 W/m^) and the detuning was set to A = —7/2. To avoid the effects of differential ground 
and excited state AC Stark shifts due to the lattice beams, the coohng and lattice beams were flashed on 
and off out of phase with each other at a frequency of 4 MHz. The rate was sufficiently large compared to 
the highest vibrational frequency encountered to avoid loss of atoms during the periods when the lattice 
beams were off. SubDoppler cooling was simulated by reducing the natural linewidth by a factor of three. 
The Doppler temperature for 6 cooling beams each with I/Ig = 0.1 and a x3 reduced linewidth is about 
54 fiK in the x, and z directions. 

In the following sections we give numerical examples of atomic compression with cooling for continuously 
variable beam angles (Sec. 13.11) and discrete angles (Sec. 13. 2p . We have used geometrical parameters that 
are compatible with future experimental demonstrations. During the compression the maximum atomic 
density is limited partly by collisions (which are not accounted for here) and primarily by radiation 
pressure from on-resonance spontaneously emitted photons. This limits the density in a magneto-optical 
trap (MOT) to a few times 10^^ atoms/cm^. However, in the lattice, the atoms are in a cylindrical pencil 
shaped trap that gets more elongated as the compression proceeds. Due to this pseudo one dimensional 
geometry the limitations due to radiation pressure are reduced. We discuss the maximum compression 
achievable given radiation pressure effects in Sec. [H 



3.1 Continuous variation of beam angles 

In this section we assume continuous variation of the beam angles using the galvanometer mirror im- 
plementation of Fig. [1^ . The mirrors are placed such that the beam separation varies from 1 mm to 
15 mm in front of a / = 50 cm focal length lens that focuses the beams in to the atomic cloud. The 
beam waists in the front focal plane are set at 250 fim in the x — y plane, and 75 /xm along z (axial 
direction) giving w^-y = 0.54 mm and Wz = 1.8 mm at the atoms. Thus in the galvanometer mirror 
setup, the half beam crossing angle can be changed from 0.057 < 9 < 0.86 deg., giving a lattice periodicity 
of 430 /im > A > 28 /xm. This is an overall compression ratio of about 15. 

The lattice beam detuning was set to A = —27v x 15 GHz and the peak intensity was 2 x 10^ W/m^ 
which results in a well depth ofUo/kB = 2 mK with an effective well depth from peak to the interwell ridge 
of 1 mK. This choice of detuning and intensity gives a scattering rate of 9.1 x 10^ photons/sec and requires 
a total power of 1.54 W for each pair of beams. This high power is a result of using both a deep well and 
our choice of beam waists. The transverse beam waist was chosen such that the atomic cloud would fit into 
a single site of the initial large period lattice. The axial beam waist was chosen so that the compressed 
cloud could be elongated along z giving a large optical depth. Relaxing any of these parameters will result 
in lower power requirements. The cooling beams give a scattering rate of 3.8 x 10^ photons/sec which is a 
factor of about 42 larger than the scattering rate due to the optical lattice, which therefore only perturbs 
the cooling weakly. 

The initial atom distribution was a Gaussian with standard deviations of (Jx = Cy = az = 250 fim. The 
initial temperature was set to = 140 /xK in x,y and z. The cloud was then placed in the potential 
defined by the first step in the lattice compression and cooled for 4 ms. The cloud was then compressed by 
continuously changing 6 between the minimum and maximum values given above over a time of 100 ms. 
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Figure 3. (Color Online) Compression with continuous change of beam angles: a)The initial Gaussian density distribution, b) the 
density after being placed in the trap and cooled for 4 ms, c) a contour plot showing the final atom distribution in the lattice, and d) 

the final number density. 
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Figure 4. The standard deviations of the atomic cloud along all three axes during the compression. 

This produced a minimum lattice site area of A(^max) — 810 jj^uP. The results of a simulation run with 
1000 atoms are shown in Fig. [3] and the change of the cloud widths with time are shown in Fig. [H The final 
compression captured 89% of the atoms in the central well with the rest of the atoms in the neighboring 
wells. The loss of atoms in the central well is mainly due to the initial placement of the Gaussian cloud in 
the optical potential. 

3.2 Discrete set of beam angles 

In this section we vary the beam angles in discrete steps as could be produced by a pixelated SLM followed 
by a lens for focusing into the region containing the atomic sample as shown in Fig. (TJ). The SLM with 
rectangular pixels is illuminated with an elliptically shaped Gaussian beam that is thinner than the height 
of each pixel along z but covers all of the SLM pixels in the x — y plane. The lenses between the SLM 
and the atom trapping region provide a Fourier transform of the field distribution so that switching a 
single SLM pixel into the transmitting state results in a beam in the trapping region with a Gaussian 
envelope along z (perpendicular to the plane of the lattice), and a sine squared envelope in the x — y 
plane. Provided that the width of the sine function is at least a few times larger than the size of the 
lattice site created at ^min the numerical results depend only weakly on the exact form of the transverse 
field distribution. Comparing the difference between the intensity pattern created by a plane wave and a 
plane wave enveloped by a sine squared function whose electric field has a FWHM = 200 fim, the only 
noteworthy difference is on the first step. For the first step, the atoms are not confined near the bottom 
of the well and the FWHM of the sine squared is smaller than the FWHM of the central well created by 
plane waves by about 48% which results in a larger number of atoms lost on the first step. After step 1, the 
percent difference of the FWHM's of the plane-wave and sine squared models, drops rapidly and is below 
1% by step 5. Since the atoms are confined near the bottom of the well after step 1, the small difference 
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Figure 5. (Color Online) Left: Density plot of the intensity of the trap along the y— axis as 6 changes by equal increments. Right: 
Atomic density at the end of a compression sequence with changing by 0.7 deg. at each step. If an atom is at a local potential 
minimum that is not the central well, it will end up far from the central well, as shown by the secondary peaks. 
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Figure 6. Left: The number of the SLM pixel turned on at each step in the compression sequence. This represents the best possible fit 
allowed by a 640 pixel SLM for a 3% change in well area at each step. Right: Percent change in well area for the compression scheme as 

a function of compression step. 



between the two potentials is negligible. We have therefore used interfering plane waves in the x — y plane 
for the numerical work, and assumed the focused beams were circular with waists of Wz — Wy — 100 fivn 
(for e.g. a beam propagating along x) in calculations of optical power requirements. 

The beam detuning was set to A = —2ti x 5 GHz and the peak intensity was 2 x 10^ W/m^ which 
results in a well depth of Uo/kB = 6 mK with an effective well depth of 3 mK. This choice of detuning 
and intensity gives a scattering rate of 7.8 x 10^ photons/sec, and requires a total power of 5 W for each 
pair of beams. This high power is a result of illuminating the entire SLM with a cylindrical beam whose 
parameters fit the width of the SLM, but only opening two pixels at a time. The scattering rate from the 
lattice beams is about 20 % of the scattering rate of the cooling beams. This ratio could be reduced further 
by using higher powers and larger detunings. 

When using a discrete set of angles there is a finite resolution with which the angle 9 can be set. If 9 
is abruptly changed by too large an amount atoms may get trapped in a local minimum of the optical 
potential, and not be swept towards the center of the lattice as illustrated in Fig. [5l In addition if the 
initial value of 9 is not small enough some of the atoms will not start in the central well and will not end 
up at the origin. To avoid these problems a large number of SLM pixels are needed so that the steps in 9 
are small. It is convenient if the amount of heating incurred at each step is approximately constant, since 
Tc can then be held constant. This will be the case if the fractional change in the area of each lattice site 
is the same for each step. In an experiment with a finite number of SLM pixels this requirement cannot be 
met exactly. For the simulations shown below we assumed the SLMs shown in Fig. [T] each had 640 pixels 
and that the range in 9 was from ^min = 0.063 to ^max = 35 deg. The pixels used at each compression step 
were chosen to approximate as well as possible a 3% change in well area between steps, as shown in Fig. 
El This corresponds to using only 158 pixels on the SLM or 157 steps of compression. 

Optical cooling is again necessary since the atoms heat as they are being compressed. The cooling 
parameters are the same as in the galvanometer mirror setup except we cool to the Doppler temperature. 
The temperature of the atoms was monitored and when either the cooling time exceeded 3 ms or all three 
Cartesian components of the temperature were below 180 /xK, the next step in compression was taken. 
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Figure 7. (Color Online) Compression with discrete variation of beam angles: a) The initial Gaussian distribution, b) the final number 

density in a region of 25 x 25 lattice sites , c) the final number density plotted to a maximum of 5% to show the atoms that have 
escaped into adjacent wells, and d) a histogram with bins defined by the lattice area A(^max) to display the final compression results. 
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Figure 8. Left: Percent of atoms captured within the minimum well area A(^max) (solid line) and percentage of atoms captured in the 

central well(dashed line). Right: density per atom during the compression. 

The Doppler temperature for 6 cooling beams each with I/Ig = 0.1 is about 160 fiK in the and z 
directions. 

The initial atom distribution was a Gaussian with standard deviations of Cx — — — 50 /xm. The 
initial temperature was set to = 140 fiK in x, y and z. The cloud was then placed in the potential 
defined by the first step in the lattice compression. The simulation involves 157 steps and 6 was varied 
between 0.063 and 35 deg. (a maximum of 70 deg. between two beams). These values correspond to using 
a Jen Optik SLM-S 640/12 with 640 pixels followed by a / = 45.6 mm focusing lens. We chose ^max = 35 
deg. to correspond to an experiment where there may only be 70 deg. of access on each side of the trapping 
region. This will produce a minimum lattice site area of A(0max) — 0.55 fim^. 

The results of a simulation run with 5000 atoms are shown in Fig. [71 The cooling and compression time 
per step was varied from a maximum of Tc ^ 3 ms at the start of the compression sequence, to a minimum 
of 100 fis at an intermediate stage, and was ~ 1.2 ms per step near the end of the sequence, giving a total 
cooling and compression time of 115 ms. The cooling time per step at the beginning of the compression 
scheme is longer since the change in well area is larger. The final compression captured 77.7% of the atoms 
in the central well with 84.8% of the atoms within a 20 ^m square around the central spot. The loss of 
atoms in the central well is due to the initial placement of the Gaussian cloud in the optical potential, the 
temperature of the atoms, and the discrete nature of the SLM. The largest loss of atoms occurs in the 
first stage of compression where the percent change in well area is about 88.8%. About 9% of the atoms 



8 



Williams & Saffman 




-3 -2-1012 a 



3rH 

2 

t 



to; 



1 

-2' 



-3L« 

-3 



IIIIIIIIHI 
IIIIIHIIIl 
IIMHIIHl 
IIMMIIIM 
IHIHIMII 




-1 1 

x(j.fnn) 



2 3 



Figure 9. Comparison of two interfering traveling waves (left) and two interfering standing waves (right) with high intensity shown in 
white. Both interference patterns are created with = tt/S giving 26 = tt/A between beams. The traveling wave produces a ID 
interference pattern while the standing wave produces a 2D lattice structure. The incident beam propagation directions are shown 

superimposed behind the plots. 



are lost in the first step, as seen in Fig. [51 This initial loss could be reduced at the expense of higher laser 
power by having more SLM pixels or deeper wells. Once the cloud is safely contained in the lattice site, 
there is a steady decrease in atoms since the effective well depth of the trap is 3 mK and the temperature 
of the cloud is kept around 175 fiK, which leads to a nonzero escape rate. Even though the total number 
captured in the central well decreases, the density of the atoms, defined as the number of atoms within 
the central well divided by the volume of the well, still increases after every step as shown in Fig [H 



3.3 Standing wave lattices 

While the traveling wave lattices described above have a relatively simple experimental implementation, 
the intensity pattern defined by Eq. ([3]) only falls to 50% between wells. This effectively wastes half of 
the well depth requiring better cooling and more powerful lasers. An alternative approach relies on retro- 
reflecting the beams to create standing waves. Keeping the same geometry as the two traveling waves case, 
but including ret roreflect ions, the result of interfering all four beams comprising the two standing waves 
is the intensity pattern 



I{x,y) cos^ [kxcos{6)] cos^ [kysm{6)] . 



(4) 



As seen in Fig. [9] this results in a two-dimensional intensity pattern, as opposed to two-traveling waves 
which give a one-dimensional lattice. Two-dimensional compression can be achieved by compressing se- 
quentially along different axes as shown in Fig. [TOl In the first step, a pair of beams with bisector along x 
is separated which compresses the atoms along the y axis, while they fill wells which are distributed along 
X. Then the process is repeated for a pair of beams with bisector along y. Since, after the first step, all 
of the atoms are located in wells along the x-axis, at the beginning of the second step, all of the atoms 
should be trapped in the central trough of the ID standing wave on the x-axis. As increases, the area 
of the wells decreases and all the atoms will be compressed to a single well located at x = y = with well 
area of A = for ^max = f • In contrast to the traveling wave case the wells are separated by regions 
of zero intensity which is a more efficient usage of laser power. Experimentally, however, a retroreflection 
arrangement is more difficult to implement. 
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4 Discussion 



The final density of tlie compressed atoms is limited by both photon rescattering and atomic collisions. 
Typical high density MOTs provide densities rig, ^ 10^^ cm~^ while dark MOTs [20] and dark optical 
lattices [21] have been used to reach rig, ^ 10^^ cm~^. Densities as high as rig, ^ 10^^ cm~^ have also been 
obtained in a highly anisotropic elongated bright MOT [22]. We therefore expect that the final density 
achievable by the galvanometer mirror scheme is on the order of rig, = 10^^ cm~^. This is a conservative 
estimate for the SLM scheme because the final trap is very anisotropic which results in greatly reduced 
radiation pressure. Using this conservative limit, it is only possible to fit a few atoms into the central well. 
If the MOT initially has too many atoms they will spread to other wells during compression resulting in 
many filled wells instead of one filled central well. The SLM compression scheme may therefore be suitable 
for loading a small number of atoms into sites of an optical lattice with sub wavelength periodicity. In this 
context the compression scheme provides an effective zoom lens with a magnification ratio of about 500 
which would be difficult to achieve with a conventional design. Additional numerical simulations show that 
it is possible to direct the atoms into a preferred lattice site, instead of the central site, by introducing a 
small phase shift between the beams at each compression step. 

The galvanometer mirror setup with a relatively small ^max is useful for creating atomic ensembles 
with large optical depth. This is of interest for studying quantum state mapping between light and small 
atomic ensembles [23]. The initial spherically symmetric MOT cloud has an optical depth, for a narrow 
beam with a width that is small compared to the MOT, of d = dz crna(x = 0, y = 0, z) where ^^(x, z) 
is the atomic density and a — 3A^/(27r) is the resonant cross section for a two-level atom with transition 
wavelength A^. Assuming a MOT with Gaussian density distribution and cr^ = Cy = az = 250 fim as in 
the previous section, and = 5 x 10^ atoms we find d = 0.044, which would imply a very weak interaction 
between a probe beam and the MOT atoms. 

However, 80% of these same 5 x 10^ atoms can be rapidly compressed into a pencil shaped lattice site to 
achieve an optical depth of d > 80, which would provide a strong interaction between the probe light and 
the atomic sample. To show this we examine the regime where the atomic temperature is small compared 
to the well depth as in our numerical simulations. The atomic density distribution inside each elongated 
well will then be Gaussian with = cry = [K{e)/Ti]yJkBTa/Uo, and a, = {w,/2).Jk^fjUo. The atomic 
density distribution can then be approximated by 



(5) 
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Figure 11. Calculated parameters of a compressed pencil shaped cloud: a) Lattice well area, b) peak density, and c) optical depth for a 

probe beam with wq = 5 /im as a function of the lattice beam crossing angle 0. The calculations assumed Wz = 0.18 cm, a relative 
temperature of ksTa/UQ = 0.05 and that 80% of the = 5 x 10^ atoms initially in the MOT are captured into the lattice site. Frame 
d) gives the optical depth as a function of the probe beam waist k;o for ^ = 0.6 deg. 



We can define the optical depth as d = — ln(Pout/^in) where Pin, Pout are the power incident on and after 
the atomic cloud. Accounting for the Gaussian transverse distribution of a radially symmetric laser beam 
and diffractive spreading as it propagates along the cloud we find 



J-L/2 Jo 



dppUaip, Z) 



1 



1 + z2/4 



-2pVK(i+zV4)] 



(6) 



where p — is the beam waist and zr = tiw^/\. This expression is valid provided the beam 

intensity is everywhere small so that there is no saturation of the absorption. When the atomic density has 
a Gaussian profile as in Eq. ([5j) we can extend the limits of the z integration to ±00 and get the analytical 
expression 



3A 



Erfc 



(7) 



where Erfc(2:) = 1 — dte~^^ is the complementary error function. 

Figure [TT] shows the calculated well area, peak density, and optical depth as a function of lattice beam 
angle 9 and probe beam waist ioi N — ^ x 10^ atoms initially in the MOT. Ki 9 — 0.6 deg. the peak 
density reaches the photon rescattering limit of about 1.5 x 10^^ cm~^, but the optical depth along z is 
d — 82.5. This is a factor of about 1800 increase compared to the value in the initial MOT, and is large 
enough to allow close to complete absorption of a photon in the atomic cloud. The radial optical depth is 
less than unity which implies that the cloud is optically thin for radially directed photons so heating due to 
photon rescattering will be minimized [22]. This approach to compression is an alternative to anisotropic 
magnetic traps [24] . 



5 Summary 

In conclusion we have analyzed protocols for compressing cold atom clouds using time-dependent optical 
lattices. The combination of a dynamically changing lattice and optical molasses can effectively create 
highly localized samples with large optical depth. While this paper focuses on two-dimensional compression. 



Compression and localization of an atomic cloud in a time dependent optical lattice 



11 



three-dimensional compression is also possible by adding two more traveling waves whose bisector is along 
z. The result is a three-dimensional lattice which compresses in the same manner as the two-dimensional 
cases studied here. Our numerical simulations show that a high percentage of atoms can be compressed 
from an atomic cloud with a standard deviation of a ~ 250 /im to a lattice site with area A ^ 35 x 35 A^. 
This results in a pencil shaped optical trap with optical depth d > 80. An atomic sample confined in 
this way is well suited for experiments requiring strong interactions between light and atomic ensembles. 
Further compression of a small number of atoms is also possible down to an area A ~ A^, provided that 
the density limits set by photon scattering and collisions are not exceeded. 

This work was supported by an Advanced Opportunity Fellowship from the University of Wisconsin, 
NSF grant CCF-0523666, and the ARO-DTO . 
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